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The wake of a flexible foil undergoing pitching oscillations in a low-speed hydrodynamic 
tunnel is used to examine the effect of chord-wise foil flexibility in the dynamical features 
of flapping-based propulsion. We compare the regime transitions in the wake with respect 
to the case of a rigid foil and show that foil flexibility inhibits the symmetry breaking of 
the reverse Benard-von Karman wake reported in the literature. A momentum balance 
calculation shows the average thrust to be up to three times greater for the flexible foil 
than for the rigid foil. We explain both of these observations by analyzing the vortex 
dynamics in the very near wake. 



1. Introduction 

In nature, wings and fins are compliant structures (see e.g. Daniel fc Combes|2002 , and 
references therein). During locomotion, these flapping appendages endure deformation, 
actively or passively, which has an important role in the propulsive and maneuvering 
capabilities of swimming ( Fish|1999 ) and flying animals (Wo otton|1992 ). A vast amount 
of work studying the effect of wing compliance in flapping wing systems has been reported 
in the literature (see Shy y et a/.|2010| for a review), and it is clear that the redistribution 
of aerodynamic forces resulting from wing deformation can have a dramatic effect in 
performance. Chordwise flexibility in particular has been shown to increase the propulsive 



efficiency in various configurations of flapping-based propulsion (see e.g. Katz & Weihs 
1978[ |Heathcote fc Gursul||2007al |Thiria fc Godoy-Diana||2010| , and a few fluid physics 



mechanisms behind the beneficial effects of passive compliance for propulsive performance 
have been identified (Eldredg e fc Pisani|20 08; Spa gnolie et a/.|2010 Ramananarivo et al. 



2011). In this paper we show that, aside from its well-documented effect on the thrust 



force, chordwise flexibility plays an important role in the stability of the propulsive 
wake. It has been shown in the literature that the mean propulsive jet produced by 



flapping motion can sometimes lose its symmetry giving rise to a deflected jet (Jones 



efal] [19981 |Lewin fc Haj-Hariri||2003| |Heathcote fc Gursul| [200761 |Godoy-Diana et al 



2008, 2009; von Ellen rieder fc Pot hos 2008] |Buchholz fc Sm its 2008] |Cleaver et aZ.|[2QT0| 



Yu et al. 2012). The mechanism of deflection is triggered by the wake vortices being 



shed as dipolar structures, instead of as the regular array of counter-rotating vortices 
of the reverse Benard-von Karman pattern characteristic of propulsive wakes produced 
by flapping motion. For the case of a rigid teardrop-shaped foil undergoing pitching 
oscillations, the phase space for the transition between different regimes has been well 
established in previous work (Godoy-Diana et al. 2008, 2009; Marais 2011). We examine 
here the effect of introducing flexibility to the problem, performing a detailed comparison 
of the vortex dynamics in the near wake of the flapping foil for two cases: the rigid foil 
(our benchmark case), and a flexible foil with the same shape but made of a compliant 
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Figure 1. Left: Top and side schematic views of the high- aspect-ratio pitching foil placed in 
the hydrodynamic tunnel. Right: photo of the flexible foil and schematic diagram. 

material. In addition to quantifying the increase in the thrust force observed when using 
a flexible foil, we demonstrate that flexibility inhibits the symmetry breaking of the 
propulsive wake. 



2. Experimental setup and parameters 

The experimental set-up closely resembles the one described in (Godoy-Diana et al. 
2008) and consists of a pitching foil of span-to-chord aspect ratio of 4:1 placed in a hy- 
drodynamic tunnel (see figure [I]). The motion of the foil is driven by a stepper motor and 
has been set to a triangular wave of varying frequencies and amplitudes. The control pa- 
rameters on the experiment are the flow velocity in the tunnel (U) and the foil oscillation 
frequency (/) and pitching angle (9). The latter determines a peak-to-peak amplitude of 
oscillation of the trailing edge that we define as A in the case of the rigid foil. As will be 
discussed below, the actual displacement of the trailing edge A e ff equals A for the rigid 
foil but varies depending on the frequency for the flexible foil. We define the following 
non-dimensional parameters: the Reynolds number Re = Uc/v based on the foil chord 
c, where v is the kinematic viscosity, fixed to Re « 1035 in the present experiments. 
The dimensionless effective flapping amplitude Ajj = A e ff/D, and a Strouhal number 
Stn — fD/U, both defined using the foil width D as characteristic length scale. It has 
been shown recently ( Godoy-Diana et aL]|2008 ) that this three-parameter description is 
needed to allow for a proper characterization of the transitions in the wake of a flapping 
foil. The amplitude-based Strouhal number St a — Stjj x Ajj that is usually found in the 
literature was in the range 0.05 < St a < 1-2 in the present experiments. 

Time-resolved measurements were performed using two-dimensional particle image ve- 
locimetry (PIV) on a plane probing the flow around the foil at mid-span. A continuous 
wave laser (DPSS 2W@532nm) and a cylindrical lens were used to produced a light 
sheet of ^lmm width in the whole imaging region ( — 10. 4D to 26. 6D and ±10. 8D in the 
x and y directions, respectively). Each recording consisted of 2100 images sampled at 
500 Hz using a Phantom V9.1 camera with 1600x1200 pixels resolution. The flow was 
seeded with 20/mi polyamide particles. The PIV computation and post-processing were 
done using a LaVision system as well as Matlab and the PIVMat Toolbox ( Moisy|2007[ ) . 
The PIV calculation was performed with 16x16 pixel interrogation windows with 50% 
overlap, giving a resolution of one vector every 0.7mm. 
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Figure 2. Left: Photo of the flexible flap where the reflection of a light sheet shows the defor- 
mation of the foil profile (top) and schematic diagram defining A e ff (bottom). Right: Effective 
flapping amplitude A e // as a function of the Strouhal number StD for =5° (+), 7.5° (x), 10° 
(*) and 15° (★) (which means A D =0.7 (+), 1.1 (x), 1.4 (*) and 2.1 (★) for the rigid foil). 



We investigate two cases, using foils with the same geometry and dimensions (indi- 
cated in figure [T]), but differing in their structural material. The first one, already used 
in previous studies, is made of acrylic glass and will be referred to as the rigid foil in the 
following. The second one, the flexible foil shown in figure [T] is made of a silicone elas- 
tomer material (vinyl polysiloxane) , crosslinked with a curing agent to obtain a Young's 
modulus of E w 30kPa. 

The deformation experienced by the flexible foil while flapping can be characterized 
using the actual amplitude A e f f (peak to peak amplitude measured at the trailing edge of 
the foil), which is larger than the amplitude A that would be observed with a rigid foil for a 
given angular oscillation. The ratio A e ff/A thus increases with the flapping frequency for 
a given value of the imposed rigid amplitude A (see figure [2j where a schematic diagram 
of this effect is also shown). The latter observation holds for the regime considered in 
the present study, however one could likely find other regimes where A e ff/A < 1. This 
is related to the problem of an optimal flexibility (see e.g. |Michelin fc Llewellyn Smithj 



2009) as a function of the forcing frequency which, although beyond the scope of the 



present work, remains to be explored experimentally. 



3. Transitions in the wake 

A description using the frequency and amplitude of the oscillatory motion as inde- 
pendent parameters has been shown recently to be the optimum framework to fully 
characterize the quasi-two-dimensional (Q2D) regimes observed in the wake of a pitching 
foil ( |Godoy-Diana et aZ.|[2008| |Schnipper et a/.|2009[ ). The phase dia gram indicating the 



regime transitions for the rigid foil in a (S^d? Ad) ma P is recalled in Fig. |3ja), where a 
typical snapshot of the mean horizontal velocity field is shown for each region of the phase 
space (snapshots of a typical vorticity field for each case are shown below in Fig. (3^-f). 
The solid line marks the transition between a Benard-von Karman type (BvK) wake 
(where the mean flow presents a velocity deficit behind the foil) and the reverse BvK 
wake characteristic of propulsive regimes (where the mean flow is a jet). The broken line 
locates the symmetry breaking of the reverse BvK street that occurs when two successive 
vortices are shed too close to each other, the closeness being determined by a ratio of the 
flapping and advective timescales (see Godoy-Diana et aT| |2009). 



The effect of adding flexibility to the flapping foil modifies the structure of the flow in 
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Figure 3. Ampli tude vs. Frequency phase diagram for the transitions in the wake for (a) the 
rigid foil (see also |Godoy-Diana et a/.|2 009) and (b) the flexible foil. Inserted images correspond 
to a typical mean horizontal velocity field for each regime. Solid and dashed lines correspond 
to transitions between regimes (see text). For the flexible case, the transition that marks the 
symmetry breaking of the reverse Benard-von Karman wake is not observed, (c-f) Snapshots of 
typical vorticity fields for each regime. 



the wake and thus the transition map (see Fig. [3p). While there is no observable effect 
in the transition between the BvK and reverse BvK regimes (i.e. the solid line does not 
move appreciably between Figs. [3^ and [3)3), the symmetry breaking of the BvK street is 
inhibited for the flexible foil in the range of parameters explored. The propulsive jet that 
is deflected for the rigid foil (in the top-right region of the (5£d, Ad) parameter space 
in Fig. |3|i) stays symmetric in the case of the flexible foil. This can be observed in more 
detail in the comparison of the mean horizontal velocity profiles for a specific point of 
the parameter space shown in Fig. [4j 

A standard momentum balance using the mean velocity field permits to estimate the 
average thrust force (see |Godoy-Diana et a/.|2008| ). Fig.[4]shows the mean drag coefficient 
Cd/Cdo surfaces in the (Str>,An) map for the two foil^fj We recall that for a given 
angular amplitude of the pitching motion, at equal flapping frequency the effective 
amplitude has been observed in the present experiments to be larger for the flexible foil 
than for the rigid one. Now if we compare two points in the parameter space at equal 
effective amplitudes (i.e. equal Ad), the calculated forces for the flexible and the rigid 
foils are indeed similar (see contours in Fig. [4|. Seeking to characterize the effect of 
flexibility however, one should consider the pitching motion as the input and the passive 
deformation of the foil as an effect, which would mean comparing points of equal angular 
amplitude (e.g. the square marker in Fig. [4| where the propulsive force is significantly 
larger for the flexible foil. Comparing forces between the rigid and the flexible foils at 
equal values of angular pitching amplitude, we can see that the thrust force for the 



f As mentioned in Godoy-Diana et al. (2008), the drag-thrust transition (the solid contour 
where Cd — 0) happens after the transitio n from Bv K to reverse BvK regime. This observation 
has been confirmed by |Bohl Sz, KoochesfahanT ( 2009| who used an alternative control volume 
analysis that takes into account the streamwise velocity fluctuations and the pressure term. 
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Figure 4. Top row: Isocontours of a mean drag coefficient Cd/Cdo surface estimated using a 
momentum-balance approach. The black contour correspond s to Cd = where the estimated 

2008). Bottom 
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profiles of the mean horizontal velocity on the downstream boundary of the control volume used 
for the thrust force estimation (solid lines) and corresponding rms (dashed lines) for a typical 
case (marked with a square symbol in the top plots) where the symmetry breaking is inhibited 
by the effect of foil flexibility. Insets: Snapshots of the mean < U x > field showing the jet flow 
in the wake of the flapping foil. 



flexible foil can be up to three times larger than the values obtained with the rigid foil, a 
feature that seems to be mainly determined by the increase in effective flapping amplitude 
resulting from the foil deformation. 

The present experiments have thus evidenced that foil flexibility has an effect on both 
the thrust force produced by the foil and the symmetry properties of the wake. In what 
follows we attempt to gain a better understanding of the mechanisms responsible for 
these observations focusing on the formation of vortices in the case of the reverse BvK 
wake, where the vorticity generated at the leading edge and along the boundary layer 
on each side of the flap merges with the tip vortex at the trailing edge and is shed as 
a single structure each half period of oscillation. The interest in the early stages of the 
vortex formation process is two-fold: on the one hand, because vortices are formed in 
the vicinity of the oscillating foil, it is during these first stages that their influence or 
footprint on the production of hydrodynamic forces will be greatest. Second, it has been 
shown that the symmetry properties of the wake are driven by the dynamics of the very 
near wake, where the formation of dipolar structures is the mechanism that triggers the 
symmetry breaking of the reverse BvK vortex street. 
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Figure 5. Four successive snapshots of the vorticity field (from left to right) showing the 
temporal tracking of vortices in the wake of a rigid foil (top row) and flexible foil (bottom 
row) . The tracked structures will be labeled vortex I (tracked with the white square) and vortex 
II (black square), see text. The case shown corresponds to StD=0.3 and = 7.5°, the latter 
corresponding to Ad =1.1 for the rigid foil and Ad =1.5 for the flexible foil. 
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Figure 6. Temporal evolution of the streamwise (x) and transversal (y) positions of a vortex 
pair obtained from the tracking shown in Fig. [5] The symbols distinguish the rigid (triangles) and 
the flexible (circles) cases. The vortex in formation still attached to the foil (vortex I) is shown 
in empty symbols while the vortex already formed, further downstream in the wake (vortex 
II) is shown in full symbols. Errorbars show the uncertainty in the position of the extrema of 
vorticity that are used to track each vortex resulting from the resolution of the vorticity field. 



4. Vortex dynamics 

For each experimental configuration, we start tracking the vortex closest to the trailing 
edge at the moment of reversal in the direction of the flapping motion, which corresponds 
to the birth of the new vortex. The vortex is characterized in terms of its circulation and 
streamwise and transverse position in the wake as in Godoy-Diana et aI] ( |2QQ9 ). In figure 
[5j a time series of vorticity fields is shown for each type of foil, rigid and flexible, for the 
same control parameters: frequency (Stn =0.3) and pitching angle (0 = 7.5°). The two 
vortices closest to the flap (counter-rotating, enclosed by squares in Fig. [5| are tracked, 
giving the time series of position shown in Fig. [6j We will refer to the vortex closest to 
the flap (negative vorticity) as vortex I and the one shed previously (positive vorticity) 
as vortex II. The first point concerns the time evolution of the streamwise position of the 
vortices: we observed systematically that in the case of the flexible foil, vortices are shed 
closer to the origin (the rotation axis of the foil) in the x direction than in the case of the 
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Figure 7. Vorticity field in the near wake for (top) the rigid foil and (bottom) the flexible foil. 
The dashed lines mark the horizontal distance between the centers of two consecutive vortices 
in the rigid case, pointing out the larger distance between the two vortices in the flexible case. 
The self-advection velocity of the dipolar structure £/di P oie is shown schematically in both cases, 
larger in the rigid case where the dipole is more compact. 



rigid foil (this can be seen in the the track of vortex I in Fig. [6| . The difference between 
the two foils is less marked in the track of vortex II. The tracks of the y position (Fig. [6]) 
show the wake deflection observed for the rigid foil in this case, while in the case of the 
flexible foil the vortices are shed symmetrically. 

The difference between the x position of vortex I in the rigid and flexible cases can 
be understood considering the deformation experienced by the flexible foil during the 
vortex formation process: when the flap changes the direction of motion, because of the 
foil deformation the tip position moves closer to the origin and so does the shed vortex. 
It follows from this observation that the two consecutive vortices of opposite sign shed at 
each flapping period, are closer in the rigid case than in the flexible case, as illustrated 
in Fig. [7| Recalling that the symmetry-breaking in the wake directly depends on the 
relative self-advection velocity of the first dipolar structure in the near wake of the foil 
( Godoy-Diana et a/.|2QQ9 ), the fact that flexibility keeps the newly shed vortex (vortex I) 
from being shed close enough to vortex II to form a vortex pair can explain the observed 
inhibition of the deflected wake. 

We study now the temporal evolution of the circulation during vortex formation. We 
define the time origin as the moment when the foil changes direction, which triggers the 
shedding of a vortex into the wake and marks the beginning of the growth process of a 
new vortex at the trailing edge of the foil. The circulation grows in time and rapidly tends 
toward a maximum value that remains relatively constant even after the next change in 
direction (which restarts the formation process for the next opposite-signed vortex). This 
evolution can be seen in Fig. [8] If we compare this circulation evolution for flexible and 
rigid foils, it appears clearly that the vortex in the wake of the flexible foil forms faster 
than in the case of the rigid foil. 
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Figure 8. Temporal evolution of the vortex circulation for the same cases of Fig. [5] Errorbars 
show the difference between two methods of calculating the circulation using either the velocity 
or the vorticity fields. 



The link between vortex formation time and the force experienced by the foil can 



be qualitatively understood using the simple approach proposed by (Protas & Wesfreid 



2003) for a 2D wake vortex street, starting from the impulse formula that associates 



the force with the time derivative of the vorticity impulse (see e.g. Saffman 1992) F = 
— Hi fn(^ x ^)°^, where the fluid density is set to unity, the flow domain Q extends to 
infinity, r — [x, y] is the position vector, Co = (d x v — d y u)k the vorticity, u and v are the 
two velocity components, and k the unit vector perpendicular to the plane of motion. 
Considering a vortex street modeled by point vortices, the projection of the impulse 
formula along the x axis, gives the following approximate expression for the drag force : 
Fd — — ^ J2i TiA^, where i denotes the z-th vortex in the wake and Ayi its distance to 
the wake centerline. The hydrodynamic force is thus on the one hand intimately linked to 
the spatial arrangement of vortices in the wake, as well as to their sign and intensity. The 
other important point is that this unsteady contribution to the force coming from the 
vortices depends on the temporal evolution of the circulation. Considering the comparison 
between the rigid and flexible foils of the present experiments, the faster formation of 
each vortex in the case of the flexible foil is therefore consistent with the observed increase 
in the propulsive force. 

In order to shed some light on the physical mechanism that determines the faster 
formation of vortices in the wake of the flexible foil we examine the velocity field in the 
near wake. We compare in Fig. [9j the velocity profile \V(t)\ at the trailing edge of the 
foil for both cases (along the white line represented on the absolute velocity fields in 
Fig. [9j left), at four different instants distributed over one foil oscillation period. The 
profiles for the rigid and flexible foils are similar, but a distinctive feature appears when 
the foil changes direction of oscillation {j^T and y§T): the velocity jet near the tip of the 
foil is clearly more intense in the case of the flexible foil. The faster growth of the vortex 
circulation observed in Fig. [5] can therefore be explained acknowledging that the flexible 
foil curvature enhances the flow around the trailing edge of the foil, locally accelerating 
the fluid that will feed the growing vortex. 

We may also note that this is compatible with another point observed in the previous 
section: in the flexible case, once the newly formed vortex is released into the wake it 
accelerates and catches up its delay with respect to the case of the rigid foil, (see the 
x position of vortex II for both cases in Fig. [6J. This transient faster initial advection 
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Figure 9. Comparison of the near wake velocity field for the flexible and rigid foils. Left: 
Modulus of the velocity fields (phase averaged, at T/12) for rigid foil (top) and flexible foil 
(bottom). Right: Profiles taken along the white line marked over the contour plots on the left 
hand side (phase averaged, at 4 different times). 

velocity of the vortices can be explained by the deformation dynamics of the flexible 
structure: in the manner of a spring, after reaching its maximum deformation the foil 
structure comes back to its initial shape, contributing to accelerating the fluid in the 
near wake. 



5. Concluding remarks 

We have shown that adding flexibility plays a strong dynamical role on the wake 
produced by a flapping foil. On the one hand, the effective amplitude obtained passively 
due to the deformation of the flexible foil while flapping (Fig. |5| leads to an increase in 
the propulsive force with respect to the case of the rigid foil (Fig. [4]). On the other hand, 
the interaction of the shed vortices with the flexible structure inhibits the trigger of the 
symmetry breaking of the reverse Benard-von Karman wake (Fig. J3|)), neutralizing thus 
the deflection of the propulsive jet that has been widely reported in the literature. The 
latter result, which is a novel observation, brings yet more evidence that wing compliance 
needs to be considered as a key parameter in the design of future flapping-propelled 
vehicles: since not only it is determinant for thrust and efficiency, but also as we show here 
because of its role in dictating the vortex dynamics that governs the stability properties 
of the wake. 

We thank gratefully D. Pradal and G. Clermont for their valuable help manufacturing 
the flexible foils and J. Bico for useful discussions. We acknowledge support from the 
French National Research Agency through project ANR-08-BLAN-0099. 

REFERENCES 

Bohl, D. G. & Koochesfahani, M. M. 2009 MTV measurements of the vortical field in the 
wake of an airfoil oscillating at high reduced frequency. J. Fluid Mech. 620, 63-88. 



10 C. Marais, B. Thiria, J. E. Wesfreid and R. Godoy-Diana 

Buchholz, J. H. J. & Smits, A. J. 2008 The wake structure and thrust performance of a 

rigid low- aspect- ratio pitching panel. J. Fluid Mech. 603, 331-365. 
Cleaver, D. J., Gursul, I. & Wang, Z. J. 2010 Vortex mode bifurcation and lift force of a 

plunging airfoil at low reynolds numbers. In 48th AIAA Aerospace Sciences Meeting, 7-10 

January 2010, Orlando, Florida, p. AIAA Paper. 
Daniel, T. L. & Combes, S. A. 2002 Flexible Wings and Fins: Bending by Inertial or Fluid- 
Dynamic Forces? Integr. Comp. Biol. 42 (5), 1044-1049. 
Eldredge, J. D. & Pisani, D. 2008 Passive locomotion of a simple articulated fish-like system 

in the wake of an obstacle. J. Fluid Mech 607, 1-10. 
von Ellenrieder, K. D. & Pothos, S. 2008 PIV measurements of the asymmetric wake of 

a two dimensional heaving hydrofoil. Exp. Fluids 44 (5), 733-745. 
Fish, F.E 1999 Performance constraints on the maneuverability of flexible and rigid biological 

systems. In Proceedings of the Eleventh International Symposium on Unmanned Untethered 

Submersible Technology, pp. 394-406. 
Godoy-Diana, R., Aider, J. L. & Wesfreid, J. E. 2008 Transitions in the wake of a napping 

foil. Phys. Rev. E 77, 016308. 
Godoy-Diana, R., Marais, C, Aider, J. L. & Wesfreid, J. E. 2009 A model for the 

symmetry breaking of the reverse Benard-von Karman vortex street produced by a napping 

foil. J. Fluid Mech. 622, 23-32. 
Heathcote, S. & Gursul, I. 2007a Flexible napping airfoil propulsion at low reynolds num- 
bers. AIAA journal 45 (5), 1066-1079. 
Heathcote, S. & Gursul, I. 20076 Jet switching phenomenon for a periodically plunging 

airfoil. Phys. Fluids 19, 027104. 
Jones, K. D., Dohring, C. M. & Platzer, M. F. 1998 Experimental and computational 

investigation of the Knoller-Betz effect. AIAA J. 36 (7), 1240-1246. 
Katz, J. & Weihs, D. 1978 Hydrodynamic propulsion by large amplitude oscillation of an 

airfoil with chordwise flexibility. J. Fluid Mech. 88, 713-723. 
Lewin, G. C. & Haj-Hariri, H. 2003 Modelling thrust generation of a two-dimensional 

heaving airfoil in a viscous flow. J. Fluid Mech. 492, 339-362. 
Marais, C. 2011 Dynamique tourbillonnaire dans le sillage d'un aileron oscillant : Propulsion 

par ailes battantes biomimetiques. PhD thesis, Universite Denis Diderot, Paris, France. 
Michelin, S. h Llewellyn Smith, S. G. 2009 Resonance and propulsion performance of a 

heaving flexible wing. Phys. Fluids 21 (7), 071902. 
Moisy, F. 2007 PIVMat: A PIV post-processing and data analysis toolbox for Matlab. Version 

1.60 17-Apr-2007. http://www.fast.u-psud.fr/pivmat 
Protas, B. & Wesfreid, J. E. 2003 On the relation between the global modes and the spectra 

of drag and lift in periodic wake flows. C. R. Mecanique 331, 49-54. 
Ramananarivo, S., Godoy-Diana, R. & Thiria, B. 2011 Rather than resonance, napping 

wing flyers may play on aerodynamics to improve performance. Proceedings of the National 

Academy of Sciences (USA) 108 (15), 5964-5969. 
Saffman, P. G. 1992 Vortex Dynamics. Cambridge University Press. 

Schnipper, T., Andersen, A. Sz Bohr, T. 2009 Vortex wakes of a napping foil. J. Fluid 
Mech. 633, 411-423. 

Shyy, W., Aono, H., Chimakurthi, S.K., Trizila, P., Kang, C.-K., Cesnik, C.E.S. 
& Liu, H. 2010 Recent progress in napping wing aerodynamics and aeroelasticity. Progr. 
Aerospace Sci. 46 (7), 284-327. 

Spagnolie, S. E., Moret, L., Shelley, M. J. & Zhang, J. 2010 Surprising behaviors in 
napping locomotion with passive pitching. Phys. Fluids 22 (4), 041903. 

Thiria, B. & Godoy-Diana, R. 2010 How wing compliance drives the efficiency of self- 
propelled napping flyers. Physical Review E 82, 015303(R). 

Wootton, R.J. 1992 Functional morphology of insect wings. Annu. Rev. Entomol. 37 (1), 
113-140. 

Yu, M.L., Hu, H. & Wang, Z. J. 2012 Experimental and numerical investigations on the 
asymmetric wake vortex structures around an oscillating airfoil. In 50th AIAA Aerospace 
Sciences Meeting including the New Horizons Forum and Aerospace Exposition, 9-12 Jan- 
uary 2012, Nashville, Tennessee, pp. AIAA Paper 2012-0299. 



